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ABSTRACT 



A representation of nonstationary narrowband random processes in terms of 
nonstationary quadrature components is proposed in a form analogous to that used 
to represent wide sense stationary narrowband random processes. The represen- 
tation is then applied to a specific Ccise in which the nonstationary narrowband 
random process is generated by the product of white noise and a deterministic 
periodic signal and then is processed by a narrowband filter. This representation 
is used in the modeling of a bi-level pulsed noise jammer which is aissumed to 
be present in a communication channel. The effect of such a jammer on a direct 
sequence, binary phase shift keyed (DS-BPSK) spread spectrum communication re- 
ceiver is evaluated and characterized in terms of the error rate performance of the 
receiver. Families of performance curves are plotted to demonstrate the effect of 
various parameters, namely signal-to-noise ratio, jammer power to signal power ra- 
tio, and processing gain, on the error rate of the complete spread spectrum receiver. 
The analysis carried out differentiates between two cases, namely fast jammers and 
slow jammers. However, the analytical tools developed make it possible to consider 
either one of the two cases without resorting to qucisi-stationary arguments as has 
been done in the past. 
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I. INTRODUCTION 



Electronic Counter Measure (ECM) techniques have gained increasing atten- 
tion in military communications. It is often necessary and important to be able to 
analyze and predict the performance of communication receivers, especially when 
operating under hostile electronic interference. In order to improve the survivabil- 
ity of tactical communication systems operating in jamming enviomments, spread 
spectrum modulation techniques have been introduced, analyzed and implemented 
in many systems in order to mitigate the effects of intentional electronic inter- 
ference. In many of the analyses carried out to date, the jammer interference 
is typically assumed to be a stationary (or wide sense stationary) random process 
with well-specified statistical characterizations (Ref l], and the natural interference 
introduced in the communication channel is generally modeled as additive white 
Gaussian noise (AWGN). However, in practice, there are many communication en- 
vironments which do not fit the model just described. Consider, for example, a 
communication system which is being jammed by a continuous wave (CW) tone 
near the transmitter operating center frequency, or by a distorted retransmission 
of the transmitter’s own signal. The interference cannot be accurately modeled 
as a stationary random process in either case. Another typical jamming scenario 
involves a jamming signal that may be pulsed between various power levels at a 
rapid rate and is therefore not stationary [Ref 2]. 

The effects of nonstationary interference on the performance of a spread spec- 
trum communication system is not known in general. It is however possible in 
many cases to characterize nonstationary interferers and determine their effect on 
specific communication systems. In fact, if the statistics of the nonstationary 
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process are periodic, such a process is said to be cyclostationary (Ref. 3]. The 
statistics and power spectrum of a cyclostationary process are normally computed 
by using phase randomization or time averaging techniques [Ref. 3]. 

In view of the importance of the use of nonconventional forms of generating 
channel interference, the aim of this thesis is to analyze the effects of nonstationary 
narrowband interference on the performance of a spread spectrum communication 
receiver. More specifically, the results of the analysis are applied to the special case 
of a narrowband bi-level pulsed noise jammer, since it is an effective form of ECM 
(Ref 3]. Although certain results concerning the same type of jammer considered 
here can be obtained by other methods (Reference 1 and 2] the results derived are 
based on the application of quasi-stationaxy assumptions which oftentimes are not 
realistic and cannot always be justified in practice. 

In Chapter II of the thesis, methods analogoxis to those used in the repre- 
sentation of narrowband random processes are used to characterize a stationary 
narrowband nonstationary random process. Such a random process is shown to 
be specified in terms of its nonstationary quadrature components. The expression 
of the autocorrelation function of the nonstationary narrowband random process 
is then established, and, all the results are used to specify the characterization 
of a bi-level, narrowband pulsed noise jammer which can be obtained by ampli- 
tude modulation of white Gaussian noise which is filtered by a narrowband filter 
centered at some high frequency fo . 

In Chapter III, analysis was carried out in order to characterize the output 
of the spread spectrum receiver consisting of a despreader, a bandpass filter and 
a correlation detector, when a nonstationary narrowband process is applied at the 
input of such a receiver. This characterization is obtained in terms of the mean 
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and variance of the Gaussian random variable at the output of the receiver’s deci- 
sion circuitry. 

Chapter IV provides a detailed analysis of the performance of a spread spec- 
trum receiver. The binary phase shift keyed (BPSK) modulation scheme was se- 
lected as the method by which digital information is transmitted and the method 
of spreading the signal bandwidth was chosen to be a direct sequence spread spec- 
trum modulation approach. Noise and jamming were assumed to be sources of 
channel interference, using the models developed in Chapter 3 for their character- 
ization. The receiver performance specification was accomplished by deriving the 
error probability or the bit error rate. The error performance was plotted versus 
various factors that take into account signal power, interference power and band- 
width spreading. The effect of the various parameters on the receiver’s performance 
was analyzed aind discussed. 

From the general representation of the narrowband nonstationary random pro- 
cess, other forms of nonstationary interference can be proposed and their effect on 
communication receivers can be derived by similar approaches. Although the de- 
rived performance of the spread spectrm receiver was developed for BPSK signal 
modulation with direct sequence bandwidth spreading, the general results can be 
modified and extended to obtain performance evaluations of receivers in which 
other modulation schemes are employed. 
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II. NONSTATIONARY RANDOM PROCESSES 



A. NONSTATIONARY NARROWBAND RANDOM PROCESSES 

A random process is said to be narrowband if its power spectral density is zero 
except for a narrow frequency region (W Hz) around a high carrier frequency (/o), 
that is, 

W « 27t/o and 27 t/<, >> 1 

If X{jt) is a sample function of such a real random process with zero mean, it may 
be expressed as: [Ref. 4] 

X(t) = Xe(t) cos 27rf ot — X, (t) sin 27rfot (2-1) 

where the power spectral density of the random processes Xc (t) and X, (t) is zero 
except for a narrow frequency region (W Hz) around / = 0. The random processes 
Xc(t) and X,(t) are referred to as the quadrature components of X(t). 

Random processes can be generally classified into stationary and nonstation- 
ary processes. If the joint probability density function (p.d.f.) of a process X(t) 
taken at arbitrary times , ta . . . fn is invariant to arbitrary time shifts for any 
integer n, such a stochastic process is said to be stationary [Ref. 5]. A random 
process is nonstationary if it does not have the above properties. A weaker form 
of stationarity is the so-called wide sense stationarity (W.S.S.) involving first and 
second order moments of a random process [Ref. 6j. The representation of Eq. 2.1 
normally implies a W.S.S. random process X(t). 

We propose here to represent a nonstationary narrowband random process in 
a form analogous to that of Eq. 2.1 and to study the statistics of the associated 
quadrature components. This is accomplished by analyzing the output of a 



11 



narrowband filter when the input is an arbitrary nonstationary process. Before 
attempting this, a convenient method of specifying a narrowband filter is developed 
first. 

1. Lowpass Equivalent of a Narrowband Filter 

Any bandpass filter can be specified in terms of low-pass filter equivalents. 
Let (t) be the impulse response of a low-pass filter and (/) be the Fourier 
Transform of /ix, (t), that is, 

h,(t) ^ H,U) = (2.2) 

where HL{f) is not necessarily symmetric about / = 0 so that is not neces- 

sarily real. Then the Fourier Transform Hb (/) of a bandpass filter having impulse 
response hs (t) can be specified as follows: [Ref. 7] 

hB [t] Hb (/) =HbU- fo) + Hb (-(/ + fo)) (2.3) 

As shown in Figure 2.1, Hb (/) is realized by summing the frequency up shifted 
(by fo) and the frequency reversed and down shifted (by /„) Fourier Transform 
of the low-pass filter. The negative sign in the term ifx,(— (/ -|- /„)) is due to the 
frequency reversion of H^ (/) described above. 

We can remark that in general (t) is a complex quantity. However, since 
the focus is on generating real narrowband random processes, we will require that' 
hs {t) is a real function of t. By inverse Fourier Transformation of Eq. 2.3 and by 
expressing (t) in the following complex form, 

hL (0 = hLr (0 + jhLi (t) (2.4) 
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_ / 
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' Hsif) 


-fo 


fo 



Figure 2.1 Graphic Representation of H^(f) and Hg(f) 
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where both hx,r(t) and hn{t) are real functions of t, we can obtain the impulse 
response hs (t) of a narrowband filter as [see Appendix A] 

hB (t) = 2Re {h^ } (2.5) 

where i2e{(*)} represents the real part of the function in brackets. It is demon- 
strated in Appendix A that h^rit) must be an even function and h^iit) must be 
an odd function if hb{t) is to be real. 

2. Output of a Narrowband Filter 

In the case of a stationary random process input, the statistics of the 
output of a narrowband filter is well-known [Ref. 4]. We will investigate the more 
general czise where the input random process to a narrowband filter is real but 
nonstationaxy. Let Y{t) be the output due to a real input signal X{t) applied to 
a narrowband filter having transfer function Hb (/) as described in Section 1 [Fig. 
2.2]. Then Y{t) can be expressed as the convolution of X{i) and Hb (t) namely 

Hb (t — oc)X{a)da 

It is demonstrated in Appendix B that Y (t) can be further expressed as 

Y{t) = 2Xc(t) cos27r/ot — 2X, (t) sin27r/ot (2-7) 





where 



/ oo 

[h^r (t — O') cos27r/oa + knit — a) sin27r/oa] X{a)da (2.8) 

• OO 

/ OO 

[^L.(t — a) cos2tt foa — htrit — a)sin27r/oa] X{a)da (2.9) 

• OO 
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Bandpass Filter 



X{t) 



1 




/ 

1 




f) 




-fo 


fc f 



^nt) 



Figure 2.2 Narrowband Filtering of Real Input X(t) 



There is a certain degree of similarity in the mathematical form of Eq. 2.7 and 
Eq. 2.1. However, the properties of Xc(t) and X,{t) given by Eq. 2.8 and Eq. 2.9 
respectively must now be established. 

3. Output Correlation Function for Nonstationarv Input 

By definition, the autocorrelation function of any random pro- 

cess y(t) is the expected value of the product of Y'(ti) and Y{t 2 ) for arbitrary time 
instants ti and [Ref. 5], that is 

Ry{ti,t2) = .E{y(ti)j/(t2)} (2.10) 



In order to obtain for the output Y{t) of the filter described in Section 

2, it will be necessary to evaluate the following quantities (since Y (t) is the linear 
combination of X^{i) and X, (t)) 



RS,,t^)^E{X,[t,)X,[t^)} (2.11) 

i2.(ti,t2)=£^{X.(<0^.(<2)} (2.12) 

RcAt^,t^)=E{XSx)X,[t^)} (2.13) 

i?.c(<l,<2) =^^{X(<i)X,(t2)} (2.14) 



The four quantities above are complicated functions of ti and <2 [See App. B]. It can 
be observed that in general Rc{ti,t 2 ) ^ R,{ti,t 2 ) and Rc,{ti,t 2 ) ^ —R,c{ti,t 2 )', 
unlike the case of stationary narrowband process, both the quadrature compo- 
nents have identical autocorrelation function and the two cross-correlation func- 
tions are null if both X^ (t) and X, (<) are real [Ref. 4]. The autocorrelation function 
Ry{ti,t 2 ) of the nonstationary narrowband process Y[t) can be formulated as fol- 
lows [ See App. B]: 



Ry{tl,t2) =4 



Rc{ii , ^ 2 ) cos 27rfoti cos 2irfot2 + R, (ti ,t 2 ) sin27r/„t sin7r/ot2 



- Re. (< 1 ,< 2 ) cos27r/„ti sin27r/et2 - R.c (^i , < 2 ) sin 27r/et cos 27r/„t2 



(2.15) 
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It is clear that i2y(f i , < 2 ) is also a linear combination of (ti , < 2 ) > n ^ 2 ) 

and R,c{ti',t 2 )‘ 



B. A MODEL FOR THE NONSTATIONARY INPUT 

The results of the previous sections will be applied now to a specific general 
class of nonstationary random processes which are generated by a system to be 
described now. This is done in order to continue the development for a class of 
nonstationaxy random processes of interest and to be able to apply the results 
to some practical problems involving jamming of spread spectrum communication 
systems. Let the nonstationary input X(t) described previously be the result of 
mixing a stationary random process W[t) with a deterministic signal q{i) cis shown 
in Fig. 2.3. We further assume that the autocorrelation function of the 

stationary random process W[t) is 

= E[W{t,)W{U)\= 6{t^ -t,) (2.16) 

That is, W[t) is assumed to be a white process with unit power spectral density 
(PSD) level. Therefore 

R.{t,M)=E[X{t,)X{U)] 



= E[W{U)q{U)W{U)q[U)] 

= q{h)q{t,)E[W{t,)W{t^)] 

= 9 ^ (^1 )^(^2 “ ^ l ) = 9 ^ ( i 2)^(^2 ” ) 



(2.17) 
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since (fi , i 2 ) is zero unless =ti. We can now apply the results from Appendix 
B in order to characterize the output of a narrowband filter when its input is a 
nonstationaiy random process of the type being described in this section. It is 
demonstrated in Appendix C that the autocorrelation terms of Equations 2.11 - 
2.14 simplify significantly when X{t) = q{t)W{t)f namely 

and 

Rc.{h,t2)=R»c{tl,t2)=0 (2.19) 

Finally, the autocorrelation function Ry[ti,t 2 ) for the nonstationary process Y{t), 
namely the output of the narrowband filter, has the following form: [see App, C] 

i,f 2 ) = 2cos27r/o(t2 — fi) / — a)hi,{t 2 — a.)q,{a)da (2. 

y — oo 

where 

q,{t)=q^t) (2.21) 

Observe that Ry{ti,t 2 ) is not a function of the time difference [t 2 — ti), so that 
Y{t) is (as expected) a nonstationary random process. 

C. A SPECIFIC EXAMPLE 

In the previous sections, we have established a representation in general form 
for a narrowband nonstationary random process generated by applying the mixture 
process q{t)W (f) at the input of a narrowband filter. As a specific example, it is now 
desired to use these results to model and characterize a pulsed noise jammer that 
is processed by a narrowband filter. Assuming that we have a bi-level pulsed noise 
jammer, we can model the jamming signal Ar(f) as the result of mixing white noise 





(2.18) 
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W{t) of unit PSD level with a bi-level deterministic periodic signal q{t). Such 
a j amm ing signal present in the communication channel would cause a typical 
receiver to produce a ja mm ing component at the output of the receiving front end 
narrowband filter. This component is denoted Y{t) as shown in Figure 2.4. 

The autocorrelation function (ti ,t2) of the input process to the narrowband 
filter is given by Eq. 2.17 while the corresponding output autocorrelation i?y (ti ,<2) 
is given by Eq. 2.20. T, is the period of the periodic signal q{t) which remains at 
the level A for the time interval [o, pT^] and at the low level C for the time interval 
[pT„T,],0 < p < 1 . [Fig. 2.5] The narrowband filter Hg (t) is based on a low-pass 
equivalent which is assumed to have impulse response hi [i) = B • e~^*U{t) where 
U{t) is the unit step function, and B specifies the 3db point of the filter. The 
details of the derivation are presented in Appendix D, where it is demonstrated 
that for the special case under consideration, Eq. 2.20 becomes 



f 42 D 

Ry{ti,t2) = 2 COS 2;r/« (t2 - 






2 



+ 



g2BT, .g2BT,l<i| _g2BpT, .g2BT,(rI 



1 _ e-2BT, 



1 



where [z] is the largest integer less than or equal to z 



r = -p+ — Min(ti,t2), 



5 = -1 + — Min(ti,t2) 

■^1 

and 



Min(ti,t2) = minimum of ti and 



( 2 . 22 ) 



(2.23) 

(2.24) 
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Figure 2.4 Schematic Representation of the Bi— level Jammer Model 





Figure 2.5 Punctians q(t) and q 3 (t) 
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Although equation 2.22 is not a simple mathematical form, it is possible to 
derive the average power of the random process Y{t), by time averaging Ry[t,t). 
The details are again presented in Appendix D, where it is demonstrated that 

Ry{t,t)=B[A^p + ^{l~p)\ (2.25) 

This result can be easily verified from the zissumptions made at the beginning of 
this paragraph. Furthermore, it can be seen that for the special case under con- 
sideration, Ry{t,t) represents the jammer power at the output of the narrowband 
filter Hsif). 

It has thus been shown that a nonstationary narrowband process can be rep- 
resented by the nonstationary quadrature components in a form analogous to that 
for the representation of wide sense stationary narrowband processes. A general 
expression of the correlation function for a filtered nonstationary input is obtained. 
The results are applied to a special case in order to model a bi-level pulsed noise 
jammer. The effects on a spread spectrum communication receiver due to such a 
pulsed noise jamming signal will be investigated in the next chapter. 
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III. EFFECT OF SPREAD SPECTRUM SYSTEM 



ON NONSTATIONARY NOISE 

A. MODEL FOR THE SPREAD SPECTRUM RECEIVER 

In spread spectrum modulation, digital information is transmitted by employ- 
ing a trasmission bandwidth which is much larger than the minimum bandwidth 
required to transmit the digital information via conventional means. In the spread 
spectrum receiver, the demodulation must be accomplished, in part, by correlation, 
of the received signal with a replica of the “spreading” signal used in the transmit- 
ter in order to recover the digital information signal. If the bandwidth is spread by 
direct modulation of a data-modulated carrier with a wideband spreading signal 
or code, the technique is referred to as direct sequence (abbreviated DS) spread 
spectrum modulation. 

The simplest form of DS spread spectrum modulation is DS-BPSK spread 
spectrum transmission where a binary phase shift keyed (BPSK) modulated car- 
rier is spread by a wideband “code signal.” [Ref. 2]. The spreading operation can 
be mathematically represented as a multiplication of the carrier by a function c(^) 
which takes on the values ±1 at a digital rate iE<, periodically. (Fig. 3.1) The trans- 
mitted spread spectrum signal is received together with some type of interference 
and/or Gaussian noise. Demodulation is accomplished in part by remodulating 
with the spreading code c(t). The process is commonly referred to as despreading 
and is a critical function in all spread spectrum systems. 

A simplified model of the DS-BPSK receiver is illustrated in Fig. 3.2. The 
figure shows the received signal r(t) first processed through a front-end bandpass 
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Figure 3.1 BPSK Spreading Code c(t) 
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Figure 3.2 Simplified DS— BPSK Spread Spectrum Receiver 



filter centered at the carrier frequency /o and having a bandwidth wide enough 
for the spreading signal (contaminated by thermal noise and perhaps a jammer) 
to be allowed through with negligible distortion. The output of the front-end 
bandpass filter is then mixed with a local synchronized replica of the code spreading 
signal c(i) in order to despread the signal. The despreading operation produces a 
contaminated BPSK signal from which the digital information can be recovered by 
conventional techniques using a filter and a correlator as shown in Fig. 3.2. 

It is well-known that in the absence of any jamming, the perfomance of the 
receiver in Fig. 3.2 in terms of the bit error probability (P(>) is given by: [Ref. 2] 

where Rb is the bit energy of the received signal, No is the one-sided power spec- 
tral density level of the AWGN interference and (5(*) is ihe complementary error 
function defined by 




B. EFFECT OF NONSTATIONARY NOISE ON THE DESPREADER 
AND BIT DETECTOR ON NONSTATIONARY NOISE 

The jamming signal produced by the bi-level pulsed noise jammer may not 
necessarily be narrowband. However, after the filtering operation performed by 
the receiver’s front-end bandpass filter the jamming signal can be considered to 
be narrowband. We therefore investigate here the effect of the despreader and the 
bit detector system on the input nonstationary narrowband random signal which 
models the jammer. 

The receiver of Figure 3.2 is now assumed to have a narrowband nonstationary 
zero mean input y(t) as illustrated in Figure 3.3. The output of the bandpciss filter 
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Figure 3.3 Simplified DS— BPSK Spread Spectrum Receiver 
for Nonatationary Narrowband Input 



is denoted by z{t) while the bandpass filter impulse response is denoted hi {t) and 
can be specified in terms of a low-pass filter equivalents denoted hi (t). The output 
of the receiver due to the input y(f) is denoted Yy. It can be seen that 



/ oo 

hi{t - T)y{T)c{r)di 

■ OO 



(3.3) 



and 



r” 

Yj= z[t) cos 27rfotdt 

Jo 

which can also be expressed as 



Yj = f \ f hi{t - T]y{T)c{T)dT dt 

Jo iJ-oo 



(3.4) 



The jamming of the transmitted signal is always being considered as an action 
of adding a nonstationary Gaussian random process to the signal. As will be seen 
in the analysis to be carried out in Chapter IV, part B, the effect of the jammer on 
the receiver’s performance is completely accounted for by evaluating the variance of 
the output of the bit detector. That is, only the variance of Yj will be established 
and the other statistical properties of Yj which have no direct influence on the 
systeihs bit-error probability will be neglected. 

The variance Oj can be shown to be given by an infinite sum of terms, namely 

<^5= Z C„ r H,(f)H,(nR,-f)P,(-/)P.,{f-nR,). 

•/ — OO 



n — oo 



/: 



Hi, (u)Hi (nRi — u) [S.(/ + /o — «) + 5.(/ — /o — u)] dudf 



(3.5) 



where iE, is the period of the nonstationary jamming signal and C„ is given by: 



= 



A^p-\r [I— p) 
p{A^-C-) 



sinriTTp 

ll 



nnp 
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; n = 0 
; n 0 



It can be seen frona Eq. 3.5 that <7^ is a function of the duration Tb of the 
information bit, the duration Tc of the spreading code, chips, and T, , the jammer 
pulsing period. However, practical spread spectrum systems utilize spreading codes 
for which the chip duration is very much smaller than the information data bit 
duration [Ref. 8], That is 



r„ << Tb or Rc » Rb 

where is just T~ ^ , the chip rate, and Rb is T~ ^ , the bit rate. It is therefore 
possible to define here two classes of jammers, namely fast jammers for which the 
pulse repetition period T, is of the same order of magnitude as T^, so that the 
jammer pulsing rate is much larger than the data rate, and slow jammers for which 
Tc is of the same order of magnitude as Tb. The effect of the jammer pulsing rate 
on the variance Cj will be examined separately for both cases. 

1. Variance for the Fast Jammer 

Consider first the Fourier Transform Pc{f) of the function Pc{f) defined 
in Eq. E-2, namely 

Pc{f)= I \os2tt = G{f - fo) + G{f + fo) (3.6) 

Jo 

where 

G(/) = ^ sine 

We can see that Pc{f) is significant for / = ±/o over a restricted frequency band. 
Similarly, Pc[{f — nR^) is significant for / = ±(/o +nRq) over a similarly restricted 
frequency band. The functions Pdf) and Pc{f — nRd are diagramed in Figure 3.4 
from which it is possible to observe that if 

fo — tiRg -\- Rb << fo — Rb 
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or equivalently, if 



2Rb « nRg 



(3.7) 



essentially no spectral overlap of the frequency functions occur. This condition 
except for n = 0 is satisfied for what has been defined as a fast jammer since 
Rg » Rb under fast Jamming. 

It is possible to conclude therefore that in the case of a fast pulsed jamming 
environment, the only contribution due to the infinite sum of terms that make up 
the variance Cj is the term n = 0 (See Eq. 3.5). Therefore 




Hl (u)Hi, (-u) [S.(/ + /o - u) + S.(/ - /„ 



— u)] dudf 



(3.8) 



and substituting in Eq. 3.8 the corresponding mathematical description of the 
functions and 5c(*) as given in Eqs. E.3, G.6, G.8 and G.17, 

results in 



OoTb 



Or = 






1 + e 



-BTc 



+ - l) 



Be-BT. g- 6T, 



B^ 



BTb 
B{B^ +b'^) 



2bTbT,{B + bY 2bTbT,{B 



_ bTb (B2 - 62 ) bTb T, (B2 - 62 )2 

+ 



(3.9) 



1\ 



The details of the derivation of this result are presented in Appendix G. 

2. Variance for the Slow Jammer 

In the case of a slow pulsed jamming environment, the condition Rg » 
Rb is not satisfied, therefore a certain degree of overlapping between the spectral 
factors that make up the expression for cr ^ will occur. While the infinite sum of 
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Figure 3.4 Functions Pc(f) and Pc(f— nRq) 



Eq. 3.5 can be approximated accurately by a finite sum, the number no, of terms 
needed in the approximation will be such that 

tiqR, » 2Rb 



or equivalently 

no » ^ (3.10) 

Kq 

In this case, the result for cTj becomes 

"5 = £ C, r H,(f)H,{kR,-f)P.(-f)P,(f-kR,) 

J — oo 

k= — no 

f He {u)Hi {kRg - u) [5c if + fo - u) + 5<, {f - fo - u)] dudf 

J — OO 

(3.11) 

It is also clear that as k increases, jC^ | decreases so that the sum involved 
can perhaps be further truncated. The computation of ay in this case becomes 
quite complicated, however using munerical integration, can be obtained quite 
accurately. 
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IV. PERFORMANCE OF THE SPREAD SPECTRUM 



COMMUNICATION SYSTEM 
A. GENERAL OUTLINE 

In Chapter III, a complete derivation for , the variance of the Gaussian 
random variable that quantifies the power contribution of a bi-level, nonstationary 
narrowband jamming signal at the output of a DS-BPSK Spread Spectrum Mod- 
ulation receiver, hzis been undertaken. A diagram of the spread spectrum receiver 
is shown in Fig. 3.1. The receiver structure is optimized based on principles of 
Maximum A Posteriori (MAP) detection for binary modulation schemes. The re- 
ceived signal and noise is first processed by the front-end bandpass filter which has 
a bandwidth just wide enough to pziss the spread signal and the in-band noise. The 
narrowband signal is then processed by the despreader which consists of a multi- 
plier mixing the received signal with a synchronozed local replica of the spreading 
code signal. A second bandpass filter having a bandwidth proportional to the in- 
formation bandwidth follows the despreader, so as to further eliminate the noise 
power being presented to the remainder of the receiver. The data signal is then 
demodulated by a correlator and finally, a decision is produced every seconds 
which is translated into 0 or 1 data bits. This generates a MAP estimate of the 
digital data at the output of the receiver. 

The performance of such a receiver is usually characterized by the probability 
of error, which is, in essence, the probability of an incorrect decision made at the 
receiver. The probability of error can be further categorized into symbol error and 
bit error probability. However, for the case of Binary Phase Shift Keyed (BPSK) 
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modulation, the two error probabilities are equivalent and provide the basis for 
specifying the receiver’s performance. 

In the next section, we investigate the performance for the case in which a 
BPSK modulated signal is spread by direct sequence methods, and the signal is 
received in AWGN interference as well as the bi-level pulsed noise jammer described 
and analyzed in previous chapters. 



B. RECEIVER PERFORMANCE FOR BPSK MODULATION 



The structure of a spread spectrum correlation reciever is shown once again 
in Fig. 4.1. In BPSK modulation, the two transmitted signals are [Ref. 4] 

l~2E~ 

Si (t) = \l ^ cos 27rfot 0 < t ^ T(, and 

'' ^ (4.1) 

2E ^ ' 

S2(t) = -y -^cos27r/„t 0<t<Tb 

where is the energy of each signal and Tf, is the duration of each signal. Since 
each signal represents a single information bit, is the energy per bit and T~ ^ is 
the bit rate. 

At the output of the second bandpass filter, the signal y{t) can be regarded as 
the summation of the transmitted signal s,(t), narrowband noise Nth{t), and the 
jamming signal yj{t) convolved with the impulse response hj{t) of the bandpass 
filter. That is. 



(4.2) 



t = 1 or 2 



Hence, 



f 

= / y{i) cos fotdt 
Jo 

= f Si{t) * hi{t) cos 2irfotdt + f Nth{t)C [t] * hj[t) cos27t f„tdt 

Jo Jo 

fT, 

+ / yj{t)c{t) * hi {t) COS 2nf„tdt 
Jo 



(4.3) 



35 




36 



Figure 4.1 DS— BPSK Spread Spectrum Receiver 



or simply, 



Y = Si+N,^+Yj 



(4.4) 



where 5<, N^h and Yj represents the three integrals in Equation 4.3 respectively. 
It is clear that Y is in fact a random variable whose statistics can be evaluated 
once certain properties of Sj , Nth and Yj are established. 

We first assume that the noise contribution at the front end of the receiver is 
zero mean, white Gaussian noise with ^ as its two sided power spectral density 
level [Fig. 4.2]. The autocorrelation function for the noise, ednoted is 






(4.5) 



where 5(t) is the Dirac Delta function. Based on these assumptions it is shown in 
Appendix H that Nth is also a zero mean, Gaussian random variable with variance 
specifiad in the sequel. Due to the fact that the receiver is linear, the noise 



contribution at the output retains its Gaussian statistics. The expression for a 



th 



is from App. H given by (Eq. H.18). 



= 



NpThk 

2 % 



(4.6) 



where 



and 



k = 0.903 



^ Rc A . . 

Gc = — = processing gam. 

Rb 



It is shown that Si and S 2 are given by [Ref. 2] 



5, = 






and — 



EbTb 



(4.7) 
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Figure 4.2 Power Spectral Density of the Thermal Noise 
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Furthermore, 10 is the quantity related to the jammer and its variance <7^ is given 
by (see Eq. H.4) 

= ^J{G,) 

where J(Gc) is defined by Eq. 4.6. To further simplify the notation, let X = 
Nth + 10> so that 

( 4 - 7 ) 

since the noise and ja mm ing are zero mean uncorrelated processes, and 

Y = Si +X (4.8) 



then X is zero mean Gaussian random variable with pdf given by [Ref. 2] 



/. m = 



1 

e X 






(4.9)' 



For the BPSK modulation scheme that is being considered here, the decision 
regions of the signal space diagram can be illustrated as shown in Fig. 4.3 [Ref 4]: 
From equations 4.8 and 4.9, it can be seen that the probability of receiving the 
signal Y given that a signal 5,- was actually transmitted is: 



Pr{r = Y/Si{t) was transmitted} = fri,.{Y J Si) = fx{Y — 5.) (4-10) 



The probability of bit error rate is then given by: [Ref. 2] 



P, = Q 




or 



P, = Q 



EhTh 

oj, +< t 2 



(4.11) 
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sin 2 7Tfot 




Figure 4.3 Signal Space Diagram for BKSP Modulation 
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where Q(') is the complementary error function defined earlier in Chapter III. 
Substituting and from Equations H.4 and H.18 results in 



or 




C. ANALYSIS OF ERROR PERFORMANCE 



(4.12) 



As derived in the previous section, Eq. 4.16 yields the probability of bit error 
Pb for a DS-BPSK spread spectrum receiver under AWGN noise interference and 
narrowband bi-level fast pulsed noise jamming. Clearly P 5 is a function of the 
power spectral density level N„ of the AWGN, the average power of the jammer 
given by BC„ and also, the processing gain G<,. It is clear from Eq. 4.11 that 
the receiver bit error probability will increase if the noise and/or jammer power 
increases as the ( 5 (-) function is monotonically decreasing. 

In order to analyze the quantity Pt, in more detail, we define the signal to noise 
ratio (SNR) to be the ratio of Eb to No and the jammer power to signal power ratio 
(JSR) to be the ratio of average jammer power Pj to the average signal power P 
[Ref. 4]. Hence, 

SNR = ^ , JSR = ^ (4.17) 

where 



Pj = BCo = 2t^RoCo (4.18) 

P = EbRb (4.19) 
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so that Eq. 4.16 can be reexpressed as 






SNR 

903 , (JSR){SNR) 

SnGc 



(4.20) 



Based on Eq. 4.20, families of curves of Pb versus SNR for fixed JSR values can be 
plotted and analyzed. However, first plots of versus JSR for different values of 
Gc under various SNR conditions are shown. Logarithm scales are used for Pb due 
to the large variation in this parameter. From Figure 4.4a to d, one can observe 
that in general, under intense jamming (i.e., JSR large), the receiver probability of 
error is also large, thus implying that a significant number of delivered data bits 
are in error. Under a fixed JSR and fixed Gc values, it can be seen that if the SNR 
is large, Pb can be made significantly small. For example in Fig. 4.4b, it is shown 
that when JSR = 15dB and Gc =100; 



Pb{a.t SNR = 4) « 10- 

and 

Pfe(at SNR = 20) « 10-2° 

However, when the jamming is very strong (for JSR > 30 dB) Pb tends to a limit of 
approximately 10- ^ regardless of the value of all other factors. It can be observed 
from Eq. 4.16 that if JSR increases indefinitely, the limit of Pb equal to 0.5 is 
reached. 

From Fig. 4.4a to d, we can also conclude that for any fixed values of JSR 
and SNR, an increasing processing gain Gc implies a decreasing value of Pb. The 
phenomenon is more prominent in a high SNR environment than in a low SNR 
environment. There is no significant improvement in Pb for SNR in the range from 
0 dB to 8 dB however a large improvement in Pb is obtained when SNR exceeds 
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Iog,o Pb log,o Pb 

-60 -40 -20 0 -40 -30 -20 




(a) 




JSR in dB 

(_b_) 

Figure 4.4 Pb versus JSR for (a) Gc=lQ 
(b) Gc=100 
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log.o Pb log,o Pb 
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(C) 




Ld) 



Figure 4.4 Pb versus JSR for (c) Gc=200 
(d) Gc=400 
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12 dB. The immediate conclusion that we can draw from this is that by increasing 
the chip rate {Rc) of the spreading code, resulting in a higher processing gain, 
the perfromance of the receiver can be improved under the jamming environment 
analyzed, however this counter mezisure loses effectiveness when powerful jamming 
is present. 

To better illustrate the effect of on Pb, another family of curves is plotted. 
In Fig. 4.5a to d the variation of P5 8is a function of JSR is shown for different 
values of and SNR. It is clearly shown that all the curves converge toward Pb 
equal to 0.5 as JSR increases. From Fig. 4.5a we can see that Pb (at JSR = 20 
dB) reduces from approximately 10“ to 10“^ ® when G^ increases from 2 to 400 
under weak signal (SNR = 4 dB) conditions, while Pb reduces from approximately 
10“°'^ to 10“^^ under the same conditions when SNR = 20 dB (see Fig. 4.5d). 
At JSR > 35 dB, no significant gain can be obtained regardless of signal strength. 
As a result of this, we can conclude that by increeising the chip rate to counter the 
narrowband bi-level feist pulsed noise jammer is not effective under strong jamming 
conditions, however, a high chip rate can reduce the probability of bit error of the 
jamming is not too severe. 

It is also possible to represent the error performance curves by plotting logjg Pb 
versus S JR, the signal power to jammer power ratio under various SNR conditions. 
Fig. 4.6 a to d show how Pb varies with SJR. The curves can be interpreted 
in essentially the same way as the previously mentioned graphs. By comparing 
Fig. 4.6 a to d with the results obtained in Reference 2, we observe that the 
performance of the DS-BPSK spreak spectrum reciever degrades under a jamming 
environment in a manner similar to the case of barrage noise jamming or partial 
band noise jamming environment presented in the reference. Furthermore, the 
receiver performance is well below the boundary of the performance for the 
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(t) 



Figure 4.5 Pb versus JSR for (a) SNR=0 dB 
(b) SNR=4 dB 
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log,(,Pb log,oPb 
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Figure 4.5 Pb versus JSR for (c) SNR=16 dB 
(d) SNR=20 dB 
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Figure 4.6 Pb versus SJR for (a) Gc=10 
(b) Gc=100 
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Cd^ 

Figure 4.6 Pb versus SJR for (c) Gc=200 
(d) Gc=400 



49 



worst case pulsed noise jamming derived in Reference 2 based on a quasi-stationary 
analysis method. Except for the weak signal condition case, the receiver perfor- 
mance approaches the above mentioned boundary. In general, by increasing the 
processing gain, the performance of the receiver improves significantly. 



50 



V. CONCLUSIONS 



In this thesis it has been demonstrated that any nonstationary random process 
that is bandpass filtered can be represented in a mathematical form analogous to 
that used to represent narrowband wide sense stationary random processes. The 
quadrature components of the narrowband nonstationary process are themselves 
nonstationary low pass random processes. In general, the autocorrelation function 
of the two quadrature components are not identical and the cross-correlation func- 
tion of the quadrature components is not necessarily an odd function. However, 
when the input nonstationary random process is real, the correlation functions of 
the quadrature components of the narrowband output process have properties sim- 
ilar to those associated with narrowband wide sense stationary random processes 
The general representation of the nonstationary narrowband random process 
has been applied as a model for a specific class of electronic counter measure 
signals. Specifically a bi-level pulsed noise jammer, in which the power of the 
jammer periodically pulses between two levels is modeled using the previously 
mentioned nonstationary narrowband random process representation and its effect 
on a DS-BPSK receiver analyzed. This type of jammer was analyzed because it is 
quite effective in the sense that it produces significant performance degradations 
on the receiver. In practical cases, the jammer pulses between two levels at a rate 
faster than that of the data bit rate. The quasi-stationary analysis method for 
receiver performance is not valid. Therefore, a methodology has been developed 
that allows rather precise evaluation of receiver performance in terms of signal, 
noise, and jamming powers as well as spread spectrum processing gain. 
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The effect of the pulsed noise jammer on the spread spectrum binary phase 
shift keyed receiver was examined in detail and it is shown that a strong transmitted 
signal is more effective at mitigating the effect of the jammer than increasing the 
chip rate of the spreading code. The pulsed noise jamming can be shown to be 
quite effective under various signal power and processing gain values, as rather 
large bit error rates result when powerful jamming is present. 

The thesis develops a methodology that can be used to evaluate the effect of 
a nonstationary narrowband random process interference on communication chan- 
nels. There are two benefits from this study. First, in terms of ECM performance 
it makes it possible to evaluate the effectiveness of such jammers, and second, in 
terms of electronic counter counter measures (ECCM), the performance of jam re- 
sistant receivers can be evaluated more accurately. Although the example chosen 
in the study involves BPSK modulation spread by direct sequence methods, results 
for other forms of nonstationary jammings on different types of spread spectrum 
receivers can also be obtained from the above results with appropriate modification. 
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APPENDIX A 



REPRESENTATION OF NARROWBAND FILTERS 



Let 






(X.l) 



be the specification of a low-pass filter (LPF), and let 



hs (t) Hb (/) = HbU-L)+Hb (-(/ -I- /.)) 



(d.2) 



where /„ is assumed to be much higher than the cutoff frequency of the LPF. 
Thus, Hb if.) is the transfer function of a narrowband filter centered at / = ±/o- 
By inverse Fourier Transformation, we have: 

ba(f)= /” Fre(/)e'""'«i/ 

*/ — oo 

= /” HB{f- r 

J — OO J — oo 

If we express 

hL[t) = hLr{t) 3 hLi[t) 
since (t) may be complex, then 



ha (t) = [harit) + (-«)] cos 27r/ot + - /i£,< (t)] sin27r/„t+ 



3 {[hLi{t] + ha i[-t)] cos 2TTfot + [hLr{-t) - Hari-t)] sin27r/of} [AA) 

Since (t) will be assumed to be real, the imaginary part in Eq. A. 4 must vanish, 
so that 

= —hn[—t) => hai is an odd function. (-^*5) 
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and 



hLr{i) = h.Lr{—t) =>• /lir IS Ell ev6n fuiiction. (^-6) 

Observe also that 



hL{-t) = hLr{-t) +JhLi{-t) = hi,r{t) ~ jhLi{t) = hl{t) (A.7) 



Therefore 



and finally 



hsit) = 2hx,r{t) cos2nfot — 2hi, (f) sin27r/ot 
hs{t) =2Re{Ai,(t)e''="*'} 



(^. 8 ) 



(>1.9) 
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APPENDIX B 



(B.2) 



OUTPUT OF A NARROWBAND FILTER 
The output y(t) of a narrowband filter due to an input x(t) is 

/•oo 

y(t) = I hsit — a)x{a)da 

J — OO 

/ OO 

[2hir{t — a) cos2tt fo{t — a) — 2hLi{t — a) sin 27r/<, (t — a)] x{a)da 

• OO 

(B.l) 

where the representation derived in App. A has been used to obtain Eq. B.l. By 
expanding cos27r/„(t — a) and sin27r/o(t — a) it is possible to obtain 

y(f) = 2Xc{t) cos 2irf„t — 2X, (t) sin27r/„t 

where 

/ OO 

[h-Lrii ^ ck) cos27r/<,Q! + hLi{t - a) sin27r f^a] X{a)da (5.3) 

• OO 

and 

■^»(0 = [ — a)cos27T foOc — hir{t — ct)sm2'Kfoa]X {a)da (5.4) 

so that a standard narrowband-type representation is possible for the output y{t). 
K X(t) is a random process that is not necessarily stationary, X<.(t) and X, (t) are 
themselves random processes with autocorrelation functions 

= I” I" 

• I R,(a,l3)dadl3 
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and 






R.{t,,t,) = E{X.(t,)X.{t,)} 

= -i f°° r Re{Ai((, -/3)e-'’"-<»+*>}i?.(a,yJ)<ia<i/J+ 
2 y — oo y — oo 

r Re[hlih-a)hr.{t,-P)e^'^’^^^^^-^^]R4a,P)dadfi 

^ J — OO •/ — oo 



(B.6) 



The cross correlation functions are given by 






= \J^ /” Irn{^hi^{t,-oc)h:^{t2-fi)e^^^^‘^^-^^^R4a,P)dadp+ 
\r /” (5.7) 

and 

^ E {X.(t,)X,(t,)) 

= J /m - q)Ai,((j J JJ,(a,^)(ia<i/3+ 

-\r r Im{hl{t, - a]h^{t^ - R,(a,p)dad/) 

(B.8) 

Observe that in general (t) and X, (t) normally labeled the quadrature compo- 
nents of the narrowband process y{t) are not stationary in general. Nevertheless, 
the derived autocorrelations and cross correlations can be used to obtain the au- 
tocorrelation of the random process of y(t), namely 
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Ry{h,t2) = E {y{ti)y{t2)} 

= Ef^2[X^{ti)cos27rfJi - X,{ti)sm27rfJi] 

•2[Xc{t2)cos2nfot2 — X,{t2) sin 2nfot2] \ 



= 4 R^{ti,t2) cos 2 tt foil cos2TTf„t2 + iZ, (fi , 42 ) sin27r/,ii s\n2'Kf„t2 + 

- Re, {t 1, t2) cos 2n foil sin27r/oi2 

- R,e{ti,t2)sm27rfoti cos2'nfot2 (-S-9) 
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APPENDIX C 



NARROWBAND PROCESS AUTOCORRELATION FUNCTION EXAMPLE 



The autocorrelation functions determined in Appendix B are now used to 
treat a specific case, namely when the input process X(t) to the narrowband filter 
is generated by the system of Fig. 2.3. The resulting output y(t) is a narrowband 
process bearing quadrature components whose autocorrelation function is given by 
(see Eq. 2.17) 



Rc{ti,t2] = ^ j Re{/ii;(ti -a)hL{t2 - } ?, (a)da+ 

1 f°° 

- j Re {hi [ti - a)hL (<2 - a)} q, {a)da 



(C.l) 



Given that (•) represents a low-pass system and q, (a) would normally vary much 
more slowly than the frequency 2/„, it can be seen that integrals involving a 2f„ 
frequency are negligible resulting in further simplifications. Thus 



Rc{ti,t2) = - Re{hl{t^ - oc)hL{t2 - a)}q,{a)da (C. 2 ) 

— OO 

and similarly 

^,(<1,^2) = ^/ - a)hL{t2 - oc)}q,{a)da = (C. 3 ) 

J — OO 

The crosscorrelations are given by 



1 f°° 

= - J Im{hL{ti, - a)hL{t 2 - a)}q,{a)da = R,^{ti,t 2 ). 



(C.4) 
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For the special case where hi{t) is real, Rc,{ti,t 2 ) and R,c[ti,t 2 ) are both equal 
to zero, so that from Eq. B.9 and Eqs. C.2 and C.3, Ry{ti,t 2 ) becomes 



Ry{ti,t2) = 4Rc{ti,t2) [cos2irfoti cos2iifot2 +sin27r/oti sin27r/„t2] 



= 4i?c (t 1 , t2 ) cos 27T/<, [t 2 - ti ) 



(C.5) 



and finally 



Ry(ti,t2) = 2 cos 27 t/o ( tj -ti) hL{ti -a)hL{t2 - a)q,{a)da 



(C.6) 
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APPENDIX D 



AUTOCORRELATION FUNCTION AND AVERAGE POWER AT THE OUT- 
PUT OF A NARROWBAND FILTER DUE TO A BI-LEVEL PULSED JAMMER 



When q{t) is a bi-level signal as shown in Fig. 2.5, it is clear that 

,.(i) = - |(A^ - c=) \U{t - kT, - fiT,) - U(t - kT, - r,)l| (r>.l) 

where U[i) is the unit step function. Let 

/ OO 

hL {ti - Oi)hL {^2 - ol)< 1, {cc)da 

• C» 

This integral must be determined in order to evaluate Ry{ti,t 2 ) as given by Eq. 
C.6. Substituting g, (t) into I[t) and utilizing the one pole filter specification for 
(t) having 3 dB cutoff B rps, we get 

— lL_llg-B(ti + t,)g2BMin(ti,t,) _ ^2^^2 _ (^2 j g- B (t , + 1, ) 



[/: 



e^^^U{t,-a)U{t^-a) ^ U{a - kT, - pT,)da ^ 2 ) 

A;= — 00 

/ OO 00 

- a)U(t2 -a) ^ U(a - kT, - T,)da] 

■ 00 

The remaining integrals in Eq. D.2 are: 



— OO 



fc=-oo ■ *- 



g2BMin(t, ,f,) _^B(k 



A: = — OO 

if (A;-|-l)r, <Min(ti,t2) 



+ l)T,j 



(D.3) 



= 0 otherwise 
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and 



oo ^Min(ti.tj) 1 Id I- ^ 

^ |'g2BMin(t.,t,) _g2B(fc+l)T,j 



fc= — oo 



fc = -oo 

if (A: + l)r, < Min(ii,t2) 



[DA) 



= 0 otherwise 



where \z\ = largest integer less than or equal to z and 



r = -e + — Min(ti,t 2 ) 



s=l + — Min(ii,t 2 ) 
■^<1 



(£>.5j 



The difference of the two integrals given by Eq. D.3 and D.4 produces 



^ i:[ 



g2B M in(ti ,t 3 ) ^2B(k-hp)T^ 



k=-oo 



1 r 

]- E [ 



g2BMin(ti,t,) _g2B(*+l)T, 



k= — oo 



(V.6j 



since 0 < p < 1, we have r > s, so that 



M 1*1 

g2BMin(ti,f,) _ ^2B M in (t i ,t , ) (D.7) 

k= — OO k= — OO 



Furthermore 



Id hi 

g2B(*+l)T, _ g2B(fe+/>)r, 

Ac= — OO k= — oo 



;2BT, .g2BT,(*l _g2BoT, .g2Br,|r] 
1 - e-2BT, 



p.8) 



Combining the above results yields (see Eq. C.6) 



Ry{ii,t2) = 2cos27r/of • J(f) 

^ (A^B 
= 2cos2wfot< — — 

B(A^ - c^) Bdi-i-t,) 



g-B(ti + t,)g2BMin(ti,t,) _ 



g 2 BMin(t„t,) r r^l 



+ 



g2Br,g2BT,[.l _ g2BpT, g2BT,|r) 

1 -e-2sr. 



(-D-9) 
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Since X[i) represents a bi-level pulsed jajnmer, the power due to the jammer 
at the output of the narrowband filter denoted by Pj can be obtained by letting 



ti = t 2 = t m the expression of i?y(ti,t 2 ) and time averaging the result. Thus 



= A‘B- B(A^ - C^)(|r] - |.|) - 






_ 2B/>T, -2B(t-T,[r|) 



(£>.10) 



which must now be time averaged in order to obtain Pj , the average jammer power. 
With the aid of Diagram D.l, D.2, and D.3 we note that all the three functions 
{[r] ~ [ 5 ]} , {t — T, [s]} and are functions of t, and. their time average 

values can be calculated as follows; 



(H - [‘^D = lim 

n — ♦ oc 

= Urn 



(1 -p)n] = 1 - p 
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(£>. 11 ) 



g-2B(t-T,(Sl) 



1 n 

= lim -— / 

n-00 nTg Jq 



% e~‘2BT, 

ue-2B(‘+‘.) = £ 



2BT„ 



(l_e-2BT,^ (£>.12) 



g-2B (t-T, (r|) — ]{Tn 



n^oo (n -f- p)T^ 

_ e-2B/>r,(i_g-2Br,) 

“ 2BT„ 



/ e~^^*dt + n / 

Jo Jo 



(D.13) 



Therefore, the final expression for Pj is 



= B[A‘i, + cr‘{i-p)\ 






-2BT, 



2BTa 



2BT„ 



(£>.14) 
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3 - 

2 - 

1 -- 



I e 

(IHtTq) (Z-^pTq) 



■> 

t 




t H-[^] 




? 

(1+prg) 



2Tg 



(b+pT,) 



t 



Figure D.l Function [r] — [s] 
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Figure D.2 Function (t— Tq[r]) 
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Figure D.3 Function (t— Tq[s]) 
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APPENDIX E 



VARIANCE OF Yj 



Since y{t) is Gaussian and zero mean random process, F/ is a Gaussian, zero 
mean random variable due to the fact that linear operations only are applied to 
y{t). [Ref 2 ] The variance of Yj, denoted by can therefore be expressed as 



follows: 



f z{t)z{v) cos 2n /of cos 27 t / oudtdu j> 

= I I Rg{t,v) 0082-^/01 cos2irfovdtdv 
Jo Jo 



(E.l) 



where 



/ oo roo 

/ hi{t - T)hi{v - -f)Ry{T,-^)Rc{T - -f)dTd-i {E.2) 

'CO •/ — oo 



and Rc{t — 7 ) is defined as the autocorrelation function of the spreading code. For 
DS-BPSK spreading modulation, Re{S) is well known, namely 



Rc{S) = 



1-^ ;|i|<n 

0 ; |i| > T. 



where T~^ is the digital rate of the spreading code and is called the chip rate. 
The power spectrum of the spreading code Sc{f) is the Fourier transform of (5), 
namely [Ref. 2 ]. 

5,(/)=r,sinc"(/re) {E.Z) 



where 



sinc(i) = 



sin 7 TZ 
irx 
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Figure E.l Autocorrelation Function and Power Spectral 
Density of DS — BPSK Spreading Code 
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A diagram of i2c(5) and Sc{f) is shown in Fig. E.l. 
A new signal (f) is defined as follows: 



i=.W = {o 



COS 2'Kfot 0 < t < Th 
; otherwise 



So that equation E.l can be rewritten as: 



/ oo 

/ 

■ oo — oo 



(t, v)P^ (t)P<. {y)dtdv 



{EA) 



It is possible to express (t, u) as the double inverse Fourier transform of 5z (/, i/), 
that is 

/ oo roo 
• oo •/ — oo 

where 

/ oo r OO 
■ oo •/ — oo 

Let Pc{f) be the Fourier Transform of Pc(0» becomes 



Oj=f j Sz{f,i') j f Pc{v)e dtdvdfdi/ 

J — CO J — oo J — oo J — oo 

/ oo ^oo 

/ S^{f,u)P.{-f)p,{i>)dfdu 

• oo •/ — oo 

Using equation E.2 to evaluate equation E.6 results in 

Sz{f,i^)= r r P,(r,^)P4r-7) /“ T /ir (t - 

J — OO J — oo J — oo J — oo 

hi (u — 7 ) 6 '^^"*'” dtdvdrd'^- 

/ OO roo 

/ Ry {‘T, l)Rc {r - 'i)dTd'-^ 

• 00 J — 00 



since 



f hj{t- T)e~^^”^^dt = 

J — 00 



(E.i) 



(E.6) 



(E.l) 



(E.8) 



(E.9) 
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and Hi{fo) is the Fourier Transform of hj[t). Therefore 



5 = r r 

J — oo J — oo 



/ OO 

• OO j — OO 



(E.IO) 



Ry{T,'l]Rc{T - l)e 



j2n f T y 2 7T 1/ TT 



drdqdf du 



The double integration involving the variables r and 7 in equation E.IO can be eval- 
uated for the special case where y[t) is the nonstationary random process described 
previously. 

In order to evaluate <7j, the double integral D{r,q) defined as must first be 
evaluated, where 



/ OO poo 

/ Ry{T,l)Rc{T,l)i 

■ 00 — OO 



j2ir/r y27Ti/7 



dTdq 



(E.ll) 



Recall from equation 2.20 that 



Ry 7 ) = 2 cos 27 t/o (r - ^)n{r, 7 ) 



[E.\2) 



where 



/ OO 

ht [t - (7 - {ct)da 

• OO 



(E.13) 



The double Fourier transform pair 7) and H{u,w) can hence be expressed as: 



/ OO poo 

/ r/(r,7)e"^''’'^"e^'*"''"'drd7 (E.14) 

*00 J — OO 



so that 



/ OO p 00 
■ OO J — OO 



' dudw 



(£.15) 
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Therefore, 



Z)(r, 7 ) = r r [ f dudw 

•/ — <X> •/ — oo L — oo •/ — oo 



Rc{t - 7 ) cos 27 t/o (r - 7 )e 



-} 2 irfT ] 2 ni/’i 



drd^ 



/ oo TOO roo roo 

/ H{u,w) / 2i24r - 7 )cos 27 r/o(r - 7 )- 

• 00 «' — OO •/— oo«/ — 00 



drd'idudm 



= r r H(u,w) r [s,{f-u+M+s,(f-u-M\ 

•/ — 00 «/ — 00 J — 00 



d'^dudw 



/ oo roo 

/ H{u,w) [Sc{f — u + fo)+ Sc{f -u - fo)]S{f - u + w ~.u)dudw 

• OO J — 00 

(^.16) 

Carrying out the last integration with respect to w, results in 

= f H{u,u + u - f)[S^{f - u + fo) + S^{f - u - fo)]du {E.17) 

J — OO 

Evaluation of H{u,w) which is defined in Equation E.14, can be accomplished 
by performing a double Fourier transformation on » 7 (r, 7 ). The single Fourier trans- 
form of 7 (r, 7 ) is 

/ OO ^00 ^ OO 

77(r,7)e"^^"“’’dr = / /il (7 ~ a) / (r - (a)da (E.18) 

- OO •/ — ’ 00 •/ — 00 

Since 



, — j 2 ir.u a 



J — OO 

transforming again the left hand side of equation E.18, results in 

/ OO TOO TOO 

/ 7(r,7)e-^2^“^e^^"”"'drd7 = / (u)g, 

' OO •/ — OO •/ — OO 

/ oo 

(7 — d'^dcx 

’ OO 



or 



if(u,rx;) =Fi(u)if^(-u;) H da (E.19) 

J — OO 



70 



However, being a periodic deterministic signal, q, (a) can be expressed in terms of 
its exponential Fourier series expasion as follows 



q.{a)= Y, 



so that 



— oo •/ — oo 



[E.20] 



Evaluation of the C„ coefficients for a particular case of interest is worked out in 
Appendix F. Combining Equations E.19 and E.20 yields 

OO 

H{u,w) = Y, CnHL{u)HL{—w)6{nRq—u + w) [E.21] 

n= — oo 

From equation E-21, it is clear that 

oo 

H{u,u + u- f) = HL{u)HL{f -u-u) Y 5 (ni2, +:/-/) (F;.22) 

n= — OO 

Equation E.17 can hence be further simplified to yield: 



/ OO 

n — — OO 

(5c (/ + fo - u) + S„{f - fo - u]]du 

OO A Q0 

= Y C„H^R,+u-f)l H^{u)H:^{f-u-u]- 

J — CO 

n= — oo 

[Seif + fo - u) + S^{f - fo - u)\du 



(F?.23) 



Observing that in the integrand of equation E.23, u is the variable of integration, 
we can define 

S‘Af,<')= r H^(<‘)H^(f-iy-u)\S,(f + U-u) + S,(f-k-u)\du (E.24) 

J — OO 
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Substituting now D{t,'^) in equation E.l, with the use of equations E.23 and E.24 
yields: 



<^]=r r CJ(nR,+u-f)HUf,l^)di,df 

J-oo J-OO n=-oo 

{E.25) 

Integrating Eq. E.25 over the variable i/, we have: 

<T?=r Yi C,H,(f)H,(nR,-f)P.(-f)P,U-nR,)Hl(fJ-nR,)df 

•/ — oo „ _ „ 

Tl— — OO 

{E.26) 

From equation E.24, H[ (/, / — nR^) can be seen to become: 

/ OO 

Hl [u)Hl {nRg -u)[5c (/+/o ~u) +S^ (/-/o -u)]du (E.27) 

■ OO 

so that finally 

<’?= E C, r H,[f)H,{nR,-f)P,{-f)P4f-nR,)- 

„ _ „ •/ — OO 

n = — OO 

/ OO 

Hl [u)Hl [nRq - u)[5<- {f + fo ~ u) + 5<,(/ ~ fo ~ u)]dudf 

{E.28) 
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APPENDIX F 



THE EXPONENTIAL FOURIER SERIES 
OF q, (t) AND AVERAGE POWER OF THE PULSED JAMMER. 

Given q, (t) as defined in Chapter II, part C and re-illustrated as shown in 
Figure F.l since q^ (t) is a periodic signal, it is possible to expand q, (t) in terms of 
an exponential Fourier series with coefficients C „ , where: 

‘PT, 



1 f fl 

dt + 

^<J L^o JpT, 



{F.l) 



Let ti = t— in the first integral and let t 2 = t — in the second integral. 



so that 

C, 



- Jl 

~ T, J.,i^ 



1 fi^-P)Tj 2 r,(i + p) 

^2gj2»n(ti + p-J-)/T,^^^ ^ _L / c^gy2^n(t,+ -«-^ 

J-{1- 



{l-p)T,/2 



2 

(F.2) 



Using properties of even and odd functions it’ is possible to simplify Eq. F.2 to 
yield 

2nnt, 



C„ = 



1 , f‘*~ 

^9 L J o 



2 Th'Kif\ - ^ .'Ow 

COS — dti + 



n • Sin 2n7rt , j 'J, L 

■ ■ 

2n7r/r, 



^9 

pr« 

’ 1 
0 -^4 



n(l-\-p )/2 f 

Jo 



COS 



dt 



V 3n^(i + p) Sin 2n'jrt2 

2C 3 ■ ' 

2ml/ Tq 






sinriTrp sinriTrfl — /?) 

= 

rniyl — p) 



riTTp 



From equation F.3, we observe that if n = 0, then 



(F.3) 



Cq — p + (l — p) 



(F4) 
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Figure F.l Function qg(t) 



In the case where n ^ 0, C„ can be reformulated as 



= e 



jnnp 



SmUTTp 
A^p 

UTTp 






sinn7r(l — p) 
n7r(l - p) 






(f.5) 



However, expanding sinn7r(l — p) and taking into account that sinnTr = 0 and 
cos riTT = (—I)", the second term of equation F.5 becomes —p^ 
we finally obtain 



-C^) 



sin mrp 
mtp 



n ^ 0 



(f.6) 



Summarizing the results above, we have 



= 



A^p + C^l-p) 
p^(A^ _ 



if n = 0 
if n 0 



(F.7) 



The average power of the pulsed noise jammer can be obtained directly from 
the coefficients C„. Alternatively, recall the expression for from equation 

C.l and let = t 2 = t so that 



Ry (t, t) = Ry (t) = 2 cos 27 t/, 



o(i-o r 

j — < 



[t - a) hi, {t - a)q, (a) do: 



or 



R 



y 




hi {t - a)q,{a)da 



{F.8) 
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Assuming specifically that = Be ^^U[t),Ry{t) can be obtained cis follows, 

Ry{t)=2 f - a)Yq,{a)da 

J — oo 



f 



= 2B^ I 



= 2B^ ^ 



^ C'^g/ 2 ;rn(t-x)/r,^^ wherex = t-a 

= — OO 

g(2B + j2n7r /T^)x 



n = — oo 
oo 



= 2B^ ^ C„ 



-{2B.+ j27m/T,)_ 

gj2n)rt/T, 



(f.9) 



2J5 + j2n%fTg 



= 2B^ 



= 2B^ 



Cq ^ ^ ^j2nnt /T^ ^ ^ ^j2nnt/t^ 






+ E^" 



2B ^ ~"2B - j2n'KlT, ^ "2B+j2n7r/Ty 

a=l '^n=l 

oo f 

Y,Re\< 

n= 1 ^ 



r'^ ^ f -]2nirt/T, 

-^ + 2 V i2e i C„‘ — : — 

2B ^ 1 2B + j2nnT, 



Since the C„ coefficients have been shown to be of the form 



C. U(A,C,P) 



{F.IO) 



Equation F.9 becomes 

^W = 2 B ’ ^ + ^^r4ua,c,p) 

n= 1 ^ 

oo 

^ + 2^/„(A,C,^)Be 



^j2nnt/T^ 



= 2B^ 



2B +j2n7rfTg 

gy(n Jrp 4 - 2 n JTt/r, ) 



n = 1 



= 2B^ 



r> “ cos ( + mrp — tan 

^+2Y,MA,C,p)^^ 



^(25)2 + (2n7T/r,)2eJ nn/BT, 

1 nff \ 

"<« ) 



(f.io) 

It is obvious that iEy (t) is periodic in f with period T, , since the average of a cosine 
signal is zero, we therefore have 






or 



Pj =R,(t)=B[A^P + C‘(l-p)\ 



(F.ll) 
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APPENDIX G 



VARIANCE (j2 FOR FAST JAMMER 



In the case of fast jamming, the variance < 7 ^ of the jammer as demonstrated 
in Eq. 3.18 simplifies to the following expression 

<75 =c„ r iir,(/)np.(/)i^ r s',(f - u)dudf (G.i) 

J — 00 J — 00 

where 

S'(x) = Sc(x + fo) + Sc{x- fo) 

In order to simplify the notation, denote the Fourier Transform and the inverse 
Fourier Transform operation by F{(-)} and respectively and let * denote 

the convolution operation. Define 






(G.2) 



where 



V{f)= r |iTi(u)ps;(/-u)<iu 

J — 00 

so that can be reformulated as follows 

=Co r \H,(f)\'‘\p,{f)\'‘vu)df 

J — 00 

Co r \Hj{f)\^\P,{f)\^ r v{t)e-^^”^^dtdf 

J — 00 J — 00 

= C„ r v{t) r \H,{f)\^\P.{fWe-‘^->'dt df 

J — 00 J — 00 

= C„ r v(t)F-‘{|/f,(/)l"|P.(/)P} dt 

J — 00 



(G.3) 



(G.4) 
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Recall from Chap. Ill, Eq. 3.17 that: 



/’.(/) = G(/-/o) + G(/ + /o) 

G(/) = ysinc 



where for reasonable values of /o cind Tf, 



\PAfW = P,U)PA-f) « |G(/ - /o)l" + |G(/ + /o)P 



(G.5) 



or equivalently 



= (§■)' - /o)n + (y)' + MP'’ (G ®) 

The inverse Fourier Transform of this function can be easily obtained from tables 
(see Ref. 9), namely 

{i-p^/ir} = T^(^) (G.7) 



where 



(i)= ; 



.-£1 



if \t\ < T, 



if \t\ > T, 

The same considerations apply for \Hi{f)\'^ since Hj (/) is the bandpass filter trans- 
fer function specified in terms of shifted lowpass equivalents. Since (/) is signif- 
icant for / = ±/o over a restricted frequency band, 



\HiUW « - M? + - /o)p 

where the lowpass equivalent Hu,[f) is given by 

b 



HUf) = 



b + j27t/ 



so that 



\HLb{f)\ 2 + ( 2 ; t /) 2 ) 



(G.8) 



(G.9) 



(G.IO) 
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Consequently, 



F~^ 4.g-y2,r/otj 



(G.ll) 



so that evaluation of F ^ {|^/ (/)|^ |Pc(/)|^ } can be accomplished by using the 
convolution theorem [Ref. 9], that is 

{\HrifW\PcifW} =F-^{\H,{f)\^}*F-^ {\PdfW} 

- j (gj2,r/or g-y2^/or^ g-6|t-r| 2,r /« ( « - r ) y2,r/o (t - r) ^ 

= ^/ cos27T/o(t - r)dr 

. -?/:■ 



(G.12) 



where 



g[t,T) = cos27r/orcos27r/o(t — t) 

1 . (C7.13) 

= - [cos 27t/o t + cos 27t /o t cos 47t/o r + sin 27t /o t sin 47t /o t] 

Zi • 

According to the value of \t — r|, Eq. G.12 can be further expressed as follows: 






(G.14) 



where: 

/. = e-"<'-’'li(jr^g(t,T)dT (G.15) 

<((t,r)i/r (G- 16 ) 

We can now proceed to the specification of v{t) which is obtained from 

v(t)=F-'U’° \Hi.(u)\-‘-S:{f-u)du'^ 

= F-‘{|ffi(/)|’}-f-‘{S,(/)} (G.17) 



= Be ( — J cos27r/ot 
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since 



F~ if) f } = Be- ^ l‘l cos 27rfJ 






Substituting Eq. G.4 with Eq. G.14 and G.17 yields 

o-j = CqB j cos2irfot{Ii + l 2 )dt 

The limit of integration can be reduced to the interval [— TcTc] since the function 
A7r/ot(/i + l 2 )dt The limit of integration can be reduced to the interval [—Tc,Tc] 
since the function A except for |t| < hence 



=BCo 



j A cos27r/ot(Ji + (G.19) 



According to the value of t, two cases can be considered as follows: 
Case (1) Tb > t > 0, 

In this case, we have 



h+l2 = 



bT. 



(G.20) 

Since Eq. G.13 shows that g{t,r) is a function of cos47r/or and sin47r/or, it can be 
seen that the integration of functions of the form e"*’’ cos47r/or and c’”’’ sin4^/or, 
as well as forms cos47r/or and fe”*’' .sin47r/ot will yield terms involving the 
factors (m^ + 167 t^/q ) and (m^ + in the denominator of the result while 

the numerators remain bounded. The numerators are proportional to m, which in 
this case equals b ox B. Since in a practical system, 

f » b and /o » B. 
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the contribution of the corresponding terms are therefore not significant, so that, 
Eq G.20 can be simplified as follows: 



hT X 

h + h ^ • -^cos 27 t /( 






1 + ;;r dr+ 



•-/'-••(-i)'--r--(-i)'’ 



(G.21) 



Eq. G.21 can be evaluated by using standard mathematics tables [Ref. 10] and 
therefore 



bT 

h + h ^ cos 27r/ot 

O 



Case (2) 0 > t > — T&, 

In this case we have 



1 t 



g6t g- bT„ 

+ 



b bTt b^n 262 Te, 






(G.22) 



II + I2 



= ^[/-n (‘ i) 

(G.23) 

By considerations similar to those presented in Case (l), Eq. G.23 can be approx- 
imated by 



A + ^2 « -e 



bT, 1 , 

fr°^ 



+ e^ 



Tt ^ 

Tb 

bt I ^-hr 



I 



so that 



T , T n f . 

A + A « cos 27T/ot 
o 



^bt 



I t C' , 
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dr + c*’‘ f e-**' 1 


^ dr 
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> / J 




(G.24) 
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We can now proceed to the evaluation of expressed by Eq. G.19 according 
to the range of t. Since Tc <Tb, 



= 



BCpbT, 

8 






+ 



BCobT, 



8 



/ 






1 + ^ j cos^ 27r/ot 
1 — ^ j cos^ 2irfot 



„bt 



1 t 



e e 

+ 



M 



1 t 

b~W,~¥r, 



e e 

+ 



-bT 

26^ 

-bTt 



- bt I 



+ e*** dt 



262 r 



( 






(^.26) 



and using cos^ 2tt fpt = j(cos47r/ot + 1) the reasoning that gave rise to the simpli- 
fications in evaluating Eq. G.20 can be used here to neglect the terms involving 
(m^ -f IGtt^/q) or (m^ + IGtt^/q)^ in the denominatpr. Thus Eq. G.26 can be 
further simplified as follows: 



BCobTt 



+ 



16 



f 



I/:- 



1 + 






,-bTt, 



6 6T, b^T, 262Tb ^ ’ 



dt 



.-Bt 
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' 6T, 62 T 262 T '' ^ J 



dt 



(G.27) 



and using standard mathematics tables (Ref. 10) results in 
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Bp Cp bTb f 1 
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1 — e 
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52 T. 



1 — e 
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52 V 
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V yJ 
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e-BTc^bT, 
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e-BTcg-bT, 



62Tb [2(5 -6)2 Te 2(5 + 6)2T 
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52+62 



B^-b^ (52 _ 62 ) 27 ; 

g-(B-6)r. g-(B+6)T^ 



■bTt 



+ 



62 Tb 



25 



2(52 + 62) 
52 -62 “ (52 - b ^)^ T , 



(G.28) 



+ 



[b-ByZ (5 + 6)2Tjj 
which by regrouping terms, obtains the final form 



<^j = 
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(■*{ 



52TbT, 5T 



1 \ 1 4 . 

1 + ^ 



+ (2e-'»2’‘ - 1) 
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(G.29) 
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APPENDIX H 



STATISTICS OF Yj AND Ntn 



We first recall the definitions of Nth a-iid Yj from Chapter IV as follows 

’'Ti 



A r” 

Yj= y{t)c{t) * hi{t) cos 2'Kfotdt {H-1) 

Jo 

Nth = f Nth{t)c{t) * hi{t)cos2'Kfotdt {H-2) 

Jo 

Observe that Yj is a zero mean random variable since the jammer and the spreading 
code c{t) are assumed uncorrelated and furthermore, c{t) has zero mean. The 
variance of Yj has been developed in Chap. 3 and can be expressed as (see Eq. 
3.19) 






2 C'oTb J ^ 

= -g— ii + (i-c 



-BT, 



>( 



1 + e 



-BT. 



+ - 1 ) 
+ 



B^Tt,Tc BTc J BTb 
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_ 6Tfc (B2 - 62 ) 6Tb T, (B2 - 62 



(if.3) 






+ 






2bT,Z{B + by 2bTt,Z{B - by \ } 

Replacing now B and 6 in Eq. H.3 by 2ttRc and 2‘KRt, respectively thereby setting 

the 3 dB cutoff of the filters appropriately, with 

A 

° R, 

defined as the processing gain, a final expression for Cj is obtained as follows 



C^Tt 
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(HA) 



where 

•Ji^c) = 1 + 
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27t(G^ — 1)2 47T 



e~ °e 
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L(G, + 1)2 (G, 



U 

c-l)^Jj 



{H.5) 
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The determination of the variance of the random variable Nt^ can be ac- 
complished in a manner very similar to that used in the determination of . Ob- 
serving however that replacing in the jammer model the deterministic priodic signal 
q(t) with a constant unit amplitude function, namely q{t) = 1, setting 
equal to (see Eqs. 2.16 and 2.17) results in a simple AWGN noise model 

having PSD level N„ /2. Thus all evaluations involving the jammer can now be used 
with the stated modificiations in order to characterize the noise contributions to 
the performance of the system. Therefore, (see Eq. E.13) 



= f /iL(r - a)/ii('7 - a)da 

J — oo 



(H.5) 



and with 






then 






M in ( T,7 ) 









{H.6) 



Thus, th.e auto correlation function Rn{t,T) of the noise Nth{t) is 



= ■'Uos2»/<,(r - ^) 



(H.7) 



From the methodology of Chap. 3, Eq. E.14 becomes 



H{u,w)= j f iVo— 

J — oo J — oo ^ 



drd'^ 



{H.S) 



Evaluation of Eq. H.S can be shown to yield 
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Therefore, the variance of Nth (see Eq. E.25) 

/ oo 

H{u,u + 1/ — f)[Sc{f — u + fo) + S'c(/ — u + fo)]dudfdi/ 

' OO 






where 



ff(u,u + i/-/) = JV,|| 



2B 



+ ( 27 tu )2 

2B 



I $iu — (u + — /)• 



(ff.ll) 



52 + (27ru^)2 

Letting f = u — f„ or f = u + fo zs appropriate in the two integrals, obtain 

<= r i^.\H,(f+MP,(f+f^)\''z{f)<‘f+ r N„\H,{f-f,)PM-f„)\-^z{f)df 

J — 00 •/ — 00 

where 



Z{f) = f 

J — < 



B 2B 



2 B2 + (27tu)2 
Z{f) can be evaluated further as follows 
B 2B ^ B 



Sc (/ — u)dudf 



(H.13) 



k°"' ■ ""-w} 

= - f e-^^*^Rc{t)e-^'^”^*dt 

2 j-„ 



{H.U) 






00 

Tc rTc 

-Bt 



e cos27t ftdt 



-/ 



e ^ cos 2 tt f tdt 

i* 



The two integrals in Eq. H.14 can be evaluated separately and with BTc = 2jr, we 
have 



m = 



where 

z'(/) 



[l + (/T.pp 



^e-“" cos2ir/r, - 1^ ((/T.f “ l) - 2/T. sin2x/r.^ 
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From Figure H.l, it can be observed that if Rf, « Rc, then Z'{f) « 1 for / within 
the interval [— In fact, Z{f) is significant only for |/| << Rf,. Hence, Eq. 
H.12 can be simplified as follows: 



+ + + \H,(f - L)PAf - LWdf 

= r* - f.)\‘df 

^ J-R, 

{H.16) 

From Eq. G.6, we have 

\PAf - /.)P « |G(/ - / - 2h)? + |G(/)|” 

where 

G(/) = (y) sincVr.. (H.lt) 

The term \G{f — 2/<,)| has significant components outside the integration region 
[— i2e,,i2b]. Furthermore from Eq. G.8 

\H,(f - /.)P « + IHlA-IW 

and again, the term |FT£,^ (/ — 2 /„)| can be dropped from further consideration due 
to reasons similar to those discussed above in the context of \G{f — 2/o)|. Thus 

< « ^ \Hl,UW y sincVnrf/ 
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Figure H.l Function Z'(f) 
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and since |i?x,fc(/)|^ has 3 dB cutoff at / = ±Rb, the integral itself can be approx- 
imated by 



Finally 


1 

/ sinc^ {fTb)df fa k = 0.903 
2Kb J-R^ 

(7?, » (H.U) 
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